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A Prediction of Brown Dwarfs in Ultracold Molecular Gas 

Bruce G. Elmcgreen^ 

ABSTRACT 

"a^^ 

0^ i A recent model for the stellar initial mass function (IMF), in which the stellar masses are 

randomly sampled down to the thermal Jeans mass from hierarchically structured pre-stellar 
clouds, predicts that regions of ultra-cold CO gas, such as those recently found in nearby 
galaxies by Allen and collaborators, should make an abundance of Brown Dwarfs with relatively 
few normal stars. This result comes from the low value of the thermal Jeans mass, which scales 
fSJ ' as Mj oc r^/pi/^ for temperature T and pressure P, considering that the hierarchical cloud 

model always gives the Salpeter IMF slope above this lower mass limit. The ultracold CO 
clouds in the inner disk of M31 have T ^ 3K and pressures that are probably lOx higher than 
J> ■ in the solar neighborhood. This gives a mass at the peak of the IMF equal to 0.01 M0, well 

0^ , below the Brown Dwarf limit of 0.08 M0. Using a functional approximation to the IMF given 

by [1 - e-[^^/^^'l'] M-'^-^^dlogM for M > Mj, which fits the local IMF for the expected value 
^ ; of Mj ^ 0.3 Mq, an IMF with Mj = 0.01 M© in M31 has 50% of the mass and 90% of the 

2r ' objects below the Brown Dwarf limit. The brightest of the Brown Dwarfs in M31 should have 

0^\ . an apparent, extinction-corrected K-band magnitude of ~ 21 mag in their pre-main sequence 

^H ' phase. For typical star-formation efficiencies of < 10%, Brown Dwarfs and any associated stars 

O ,■ up to ~ 2.5 Mq should not heat the gas noticeably, but if the IMF continues up to arbitrarily 

JL , high masses, then the star formation efhciency has to be < 10^** to avoid heating from massive 

5-H ' stars. 

> 

t^ . 1. Introduction 

The recent observation of ultracold CO in the inner disk of M31 by Allen et al. (1995), and the 
inference that ultracold CO exists in some spiral arm dustlanes where neither CO nor HI have been detected 
in spite of a large dust extinction (M83: Allen, Atherton & Tilanus 1985; Tilanus & Allen 1993; M51: 
Tilanus & Allen 1989; MlOO: Rand 1995; M81: Allen et al. 1997) lead us to wonder about the form of the 
initial mass function (IMF) for any stars that are born at such low temperatures. Here we propose that star 
formation in ultracold gas is biased towards Brown Dwarfs, and we show that these objects, along with the 
stars that are likely to form with them, could have escaped detection up to now because they do not heat 
the gas enough to raise the temperature to normal levels, nor do they appear as stellar objects in existing 
images. 

In the Solar neighborhood. Brown Dwarfs are so uncommon that the IMF cannot continue to rise 
below about 0.1 Mq (Zuckerman & Becklin 1992; Pound & Bhtz 1993, 1995; Reid & Gazis 1997a,b). The 
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number of possible Brown Dwarfs is consistent with an extension of the IMF flattening seen between ^ 1 
Mq and 0.4 Mq (Reid 1998). The relative number of Brown Dwarfs could be much higher in the ultracold 
gas of M31, however. The detection of numerous Brown Dwarfs there would be an important confirmation 
of recent IMF models. 



2. The Case for Uniformity in the IMF 

The IMF is observed to be remarkably uniform from region to region in any one galaxy (see reviews in 
Massey 1998; Elmegreen 1999a), covering star formation that spans a factor of ~ 200 in density (Massey & 
Hunter 1998; Luhman & Rieke 1998) and a factor of ^ 10 in metallicity (Freedman 1985; Massey, Johnson 
& DeGioia-Eastwood 1995). It is also about the same in different HIT regions in various galaxies (Bresolin 
& Kennicutt 1997), and in many different galaxies on average (Kennicutt, Tamblyn & Congdon 1994), as 
indicated by the equivalent widths of hydrogen emission lines. Detailed studies of color magnitude diagrams 
in the LMC and local dwarf galaxies give the same IMF too (Greggio et al. 1993; Marconi et al. 1995; 
Holtzman et al. 1997; Grillmair et al. 1998). 

There is apparently some uniformity in the IMF with time also, because stars with a wide range of ages 
in our Galaxy all have about the same function, as suggested by halo stars (Nissen et al. 1994), and globular 
cluster stars (De Marchi & Paresce 1997). Similarly, a nearly universal IMF was found from abundance 
ratios (e.g., Fe/0, reflecting the ratio of low mass to high mass supernova processing) in QSO damped Lya 
(Lu et al. 1996) and Lya forest (Wyse 1998) lines, the intracluster medium (Renzini et al. 1993; Wyse 
1997, 1998; but see Loewenstein & Mushotzky 1996), and elliptical galaxies (see review in Wyse 1998). 

An IMF biased towards high-mass stars has been suggested for starburst regions, based on the ratio 
of luminous to dynamical mass (Rieke et al. 1980, 1993; Kronberg, Biermann, & Schwab 1985; Wright et 
al. 1988), galactic evolution models (Doane & Matthews 1993), spectroscopic line ratios (Doyon, Joseph, & 
Wright 1994; Smith et al. 1995), and infrared excesses (Smith, Herter & Haynes 1998). However, a lower 
extinction correction for M82 makes the IMF there normal (Devereux 1989; Satyapal et al. 1995, 1997), and 
more recent evolutionary models (Schaerer 1996), multiwavelength spectroscopy and broad-band infrared 
photometry (Calzetti 1997), and emission line spectroscopy (Stasihska & Leitherer 1996) give normal IMFs 
too. Large IMF shifts in starburst galaxies should also produce unobserved red populations of stars after 
the turnoff age reaches the stellar lifetime at the truncation mass (Chariot et al. 1993), and too high an 
oxygen abundance (Wang & Silk 1993). 

An IMF shift towards lower mass stars has been reported for the extreme fleld by Massey et al. (1995). 
However, this result could also come from a normal IMF in each star-forming region if massive stars form 
preferentially in high mass clouds and stop further star formation when they do (Elmegreen 1999b). The 
extreme field IMF cannot be typical, because the IMFs in most clusters and associations are about the 
same as the galaxy-integrated IMFs. 

There are also IMF dips, gaps, and a ±0.5 cluster-to-cluster variation in the power law slope (Scalo 
1998), but such variations are expected statistically given the small numbers of stars that are usually 
included in cluster studies (Elmegreen 1999b). 

There are few theoretical predictions of IMF variations. Fabian (1994) predicted that the IMF would 
be biased towards low mass stars at the cores of galaxy cluster cooling flows because of the low Jeans mass 
that results from the expected high pressures there. Larson (1999) proposed that a top-heavy IMF in the 



early Universe could explain the G-dwarf problem, the high temperature and high metal abundance of 
intracluster gas, and the large luminosities of young elliptical galaxies. Neither of these predictions have 
been directly confirmed, however. 

In what follows, we discuss the possible dependence of the minimum stellar mass on cloud temperature 
and pressure. First we discuss the minimum mass on general terms, and then we apply the results to 
ultracold clouds in M31. 



3. The Thermal Jeans Mass as a Limit to Stellar Mass 

A lower mass limit is required for star formation because cloud pieces in pre-stellar clouds extend down 
to masses much smaller than the smallest stellar mass. Clumps with masses as low as 10"'' M0 have been 
observed in great abundance as part of the normal, power-law clump mass function in the Polaris spur 
(Heithausen et al. 1998; Kramer et al. 1998). This means that the stellar mass range is only a small part 
of the total mass range for cloud clumps, and that there must be some physical process which limits the 
stellar mass at the low end. 

One possibility is that the minimum stellar mass is proportional to the thermal Jeans mass, 

Mj = 0-35(-^) (j^) Mo (1) 
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for temperature T and cloud-core pressure P. In this model, cloud pieces much smaller than Mj are not 
likely to become stars because they are not strongly self-gravitating (Larson 1992; Elmegreen 1997, 1999b). 

The mass of a star that actually forms in such a minimum unstable cloud piece can be smaller than 
Mj, perhaps by a factor of 10, because not all of the gas goes into stars, and because the resulting stellar 
system could be binary, in which case the stars would have to share the clump mass. For this reason, Mj, 
or some factor of order unity times Mj, is identified with the mass at the break point in the IMF, where the 
power law first becomes flat, not with the minimum mass of a star. This is because the model predicts only 
a fundamental change in the physical properties of cloud clumps at this mass, but does not specify how the 
stars actually form inside the clumps. Indeed, stellar masses are often observed to continue down to at least 
one-quarter of the mass of the break point, but such stars are not as common as would be expected from 
an extrapolation of the Salpeter power-law function. 

The expression for Mj is the Bonner-Ebert condition for stability of a non-magnetic, pressure-bounded 
isothermal sphere. Magnetic forces increase this critical mass (e.g., Mouschovias & Spitzer 1976), so they 
should not be included in the minimum value. Also, the mass-to-flux ratio varies in a turbulent cloud as a 
result of random compressions and magnetic diffusion, so regions with large mass-to-flux ratios will develop 
spontaneously and systematically over time, suggesting again that equation (|l|), without magnetic fields, is 
the proper condition (see also Nakano 1998). 

The numerator in Mj contains the thermal temperature, rather than the turbulent rms speed, because 
turbulence increases the critical mass, and because the smallest stars will generally form on such small 
scales that the relative turbulent speed is less than the sound speed anyway. The pressure in the expression 
for Mj is the total pressure at the boundary of the isothermal sphere. This pressure is not well defined for 
a real cloud because it comes partly from turbulent interclump motions and partly from thermal pressure in 
the interclump and overlying media (e.g.. Hunter & Fleck 1982; Ballesteros-Paredes, Vazquez-Semadeni & 
Scalo 1999). The value of it for any particular cloud is not well defined either because pressure varies from 



place to place. Nevertheless, the total pressure is much more uniform in a cloud than either the turbulent or 
the thermal pressures alone, which interchange roles during transient compressions, and also more uniform 
than the density because of the Larson (1981) scaling relations, which make pAu^ ~ constant for density p 
and rms speed Aw. In Taurus-like clouds, P ~ 10^ k^, in GMCs, P ^ 10^ k^, and in cluster-forming GMC 
cores, P ^ 10^ k^. 

These differences in P suggest that the IMF may differ as well, which means that the lowest-mass stars 
might be smaller in GMCs than in Taurus-like clouds. Any such differences would not be large, however, 
because of the square root dependence of Mj on P. Small IMF differences would not have been recognized 
either. Most stars form in GMCs under similar conditions, so their local IMFs are similar. Where the 
conditions differ, as in Taurus, there are too few stars to get an accurate IMF. In addition, the expected 
decrease in the mass at the peak of the IMF for dense GMC cores is offset by stellar mass segregation (a 
concentration of massive stars near the center) and by a slightly increased temperature from nearby OB 
stars. 

Strong variations in Mj are not expected on a galactic scale either. The numerator in Mj is roughly 
proportional to the cloud cooling rate and the denominator is roughly proportional to the background 
galactic heating rate from stars and cosmic rays (Elmegreen 1997; 1999b). Thus molecular clouds in 
equilibrium will all have about the same Mj, as long as the galactic mass-to- light ratio is about the same 
(expected variations in the IMF with M/L are discussed in Elmegreen 1999c). 

The transition from the power law part of the IMF to the low-mass turnover should contain important 
information about the star formation process, unlike the power-law part itself, which theory suggests 
contains only information about hierarchical cloud structure, i.e., about the initial conditions for star 
formation. The mathematical form for this transition, i.e., the way in which the lower mass limit actually 
affects the shape of the IMF, cannot be determined from theory yet. We have approximated this form 
simply by writing a probability P/ that a cloud piece fails to form a star on a time scale comparable to its 
internal crossing time. The functional form of this probability is unknown (cf. eqn. 0), so the exact value 
of the minimum critical mass cannot be known yet either. Differences of a factor of three in the definition 
for AIj can easily be compensated by differences in the form of Pf. Nevertheless the concept that thermal 
pressure in a star-forming cloud ultimately limits the mass of a star is an essential part of the model, and 
leads to a characteristic mass at the low end of the IMF that should scale with T^/P-^^^. 

The lower mass limit comes from another source in an alternative IMF model, which follows from ideas 
by Larson (1982) and Shu, Adams, & Lisano (1987), namely, that protostellar winds limit the accretion 
onto a star and thereby set the final stellar mass at a value close to the lower limit necessary to drive such 
a wind (Nakano, Hasegawa, & Norman 1995; Adams & Fatuzzo 1996). The wind-limited mass is not the 
same as the deuterium-burning limit itself, which is about 0.018 M© (D'Antona & Mazzitelli 1994). The 
lower stellar mass in this theory comes from a combination of the deuterium limit and the accretion rate, 
a^/G, for sound speed a (Shu, Adams, & Lisano 1987). Higher accretion rates require stronger winds and 
more massive stars before the accretion stops. This sensitivity to accretion rate makes the lower mass limit 
depend on temperature in an analogous way to Mj^ so that lower temperature clouds are likely to produce 
lower mass stars in the wind-limited accretion model too. The exact form of this temperature dependence 
and whether there is a pressure dependence are not known. 

An important difference between these two models is that the random-sampling model assumes that 
stellar mass is determined mostly by the pre-stellar clump mass. Then the IMF comes simply from random 
selection of these clumps in a certain order (densest ones first). There is, in fact, growing evidence that the 



stellar mass is determined by the pre-stellar clump mass (Casoli et al. 1986; Myers, Ladd, & Fuller 1991; 
Myers & Fuller 1993; Found & Blitz 1993, 1995; Motte, Andre, & Neri 1998; Testi & Sargent 1998), and 
not by wind-limited accretion in a large reservoir. The assumption of stellar mass proportional to clump 
mass is reasonable if even a modest fraction, e.g., > 10%, of a clump's mass goes into the stars it forms, 
because then the total variability in possible stellar mass for a given clump mass is much less than the total 
mass range of the IMF, which is a factor of 1000. In this case most of the IMF has to come from clump 
mass variations and not from variations in the star-to-clump mass fraction. Moreover, even if a random and 
small fraction of a clump's mass gets into a star and the rest goes away in a wind or dispersed disk, the 
stellar mass will still be proportional to the clump mass and the IMF will be unchanged (Elmegreen 1999b) 
if the average of this fraction is about constant. 

Given these differences between the two IMF models, the primary way to distinguish between them 
is to determine the fraction of the clump mass that goes into stars. If this fraction is either constant or 
relatively large (larger than the inverse of the total mass range for stars), then the IMF must be tracing the 
pre-stellar clump spectrum in some way. 



4. A Prediction of Brown Dwarfs in Ultracold Gas 

IMF models with a minimum stellar mass that depends on the thermal temperature and pressure in 
the star-forming cloud predict that ultracold gas in the inner disk of M31 (Allen et al. 1995) should produce 
a great abundance of Brown Dwarfs. This gas is located in the inner disk of M31, in the midst of the large, 
dark dust features studied by Hodge (1980). The surface density of stars is comparable to or higher there 
than it is locally because it is at a smaller galactocentric radius than the Sun, and the surface density of 
the dark CO clouds is higher than the local average too. Thus the pressure in the cloud cores, and perhaps 
in the average interstellar medium at this radius too, is larger than it is locally, by the product of the ratios 
of these two surface densities. 

This pressure change lowers Mj by about a factor of 3 compared to the value of Mj for local star 
formation in molecular clouds. Usually such a pressure increase is accompanied by a larger stellar volume 
emissivity and cosmic ray flux in these inner disks, both of which increase T along with P. But in M31, the 
stellar blue surface brightness and the non-thermal radio continuum flux are lower than in the inner Milky 
Way, probably as a result of the earlier Hubble type for M31. This implies that radiative and cosmic ray 
heating rates are low, so the molecular gas equilibrates at a lower temperature, around 3.5 K (Allen et al. 
1995; Loinard & Allen 1998). In that case, there is a factor of ^^ 3 drop in the temperature of potential 
star-forming material for the ultracold CO gas in M31, and, when combined with the expected higher 
pressure, a resulting net decrease in Mj by a factor of ^ 30. This gives a value of 

Mj - 0.01 Mq (2) 

in the inner disk dark clouds of M31. If the pressure is the same in the cloud cores as it is locally, then 
Mj ^ 0.03 M0. Observation of a significant population of ~ 0.01 — 0.1 Mq stars and Brown Dwarfs in the 
inner dark clouds of M31 would therefore support theories of the IMF based on either the thermal Jeans 
mass limit or the thermal accretion rate limit for the lowest mass star. 

The random sampling model also predicts that the slope of the IMF above the thermal Jeans mass 
should be independent of T, P and other cloud properties since clouds presumably partition themselves in 
a standard way by turbulence (for a discussion of cloud structure and turbulence, see Falgarone, Phillips 



& Walker 1991). We expect the same Salpeter slope, namely, —1.35 on a log — log plot, for stars with 
M >> 0.01 M0 in the M31 inner dark clouds as we find in the intermediate to high mass parts of the IMF 
measured locally. A simulation showing this independence between the high-mass slope and Mj was shown 
by figure 4 in Elmegreen (1997). Thus the random sampling model predicts that the IMF in ultracold clouds 
should be exactly shifted towards lower mass without any change in the slope of the power-law part. 

The observed IMF can be approximated by the function 

niog{M)d\ogM w A (1 - e-(^^/^^^)°) A/^^'^^dlogM (3) 

for constant A and for a in the range from 1 to 2; a ^ 1 for an IMF that flattens at low mass, and a > 2 
for an IMF that turns over at low mass. There are only a few observations yet of a low mass turnover on a 
log — log plot (Reid & Gazis 1997; Hillenbrand 1997; Nota et al. 1998), but numerous observations show a 
flattening. Note that there are usually many stars in this flattened region, and that the minimum stellar 
mass can be as low as one-quarter the mass at the turnover point. In the model, a appears only in the 
probability for failure to form a star {Pf ex e"'^*^/*^-') ). 

Equation (H) can be integrated to give the total stellar mass in various mass ranges. If Mj = 0.01 M© 
and the Brown Dwarf limit is Mbd ^ 0.08 M0, then 0.479 of all the mass M > Mj will be between Mj 
and Mbd- Approximately half the stellar-like mass down to the peak in the IMF would be in the form of 
Brown Dwarfs. 

5. Observable consequences of a Brown Dwarf bias in ultracold star formation 

The formation of a sufhciently large density of Brown Dwarfs in an ultracold CO cloud could raise 
the gas temperature and make it more "normal." The energy input required to heat a molecular cloud to 
~ lOK is about 10~^^'^no erg cm~'^ s~^ for molecular density uq (Neufeld, Lepp & Melnick 1995). This 
corresponds to a ratio of luminosity to mass equal to 6.4 x 10"'* erg s~^ g~^. If the ratio of the stellar 
luminosity that heats the gas divided by the gas mass, (L/M) from Brown Dwarfs and other young stars 
is close to this, then the gas would not be supercold. A typical L/M for pre-main sequence Brown Dwarfs 
is ~ 1 in these units (D'Antona & Mazzitelli 1994), so if the star formation efficiency is 1% and a high 
fraction of the stellar radiation heats the gas, then the stellar luminosity per unit total mass would be 0.01 
erg s~^ g"^, which is 16 times larger than the ratio needed to heat the cloud to lOK. 

Most embedded stellar radiation does not go into the gas directly, though, it goes into the dust which 
radiates it away in the IR. The dust luminosity scales approximately with dust temperature as T^ (e.g., 
Hollenbach & McKee 1979), and a typical cold dust temperature in M31 is ~ 16 K (Haas et al. 1998). If 
the dust and gas are thermally coupled in ultracold CO clouds, then the dust luminosity has to be less in 
~ 3K clouds than it is locally by a factor of about 5^ ^ 3000. The background radiation field in M31 is 
brighter than 1/3000 times the local value, so the dust and gas cannot be well coupled in the M31 ultracold 
CO clouds. This means that the average gas density in these clouds is much less than ~ 10* cm"^. The 
density will be higher than this in star- forming cores, but these are shielded from outside light by dust, so 
the core dust temperature should be low there. 

The likely range for the average cloud density can be estimated from the inclination-corrected column 
density of ~ 100 Mq pc~^ (Loinard & Allen 1998). Considering a typical cloud projected size of ~ 100 pc, 
the cloud thicknesses are probably somewhere between '^ 1 pc if they are thin shells, and ~ 100 pc, the 
thickness of the galaxy. This puts the average cloud density between 1 and 100 M© pc~"^, which corresponds 



to a molecular hydrogen density of 15-1500 cm^^. This is consistent with the result of Loinard, Allen & 
Lequeux (1995), who estimated a density of ~ 100 cm~^ from CO line ratios. Thus the average density is 
indeed too low for thermal coupling between the gas and dust. 

A better way to determine if embedded Brown Dwarfs can significantly heat ultracold CO is to compare 
the summed luminosities of these stars to the incident and embedded luminosity from field stars. If the 
Brown Dwarf luminosity in the clouds is less than the total absorbed field star power, then the clouds would 
not show any excess emission in either dust or gas from the embedded young stars. 

The total field star power received by a spherical cloud of radius R is 

Pf = AirR^jf (ttA + i?/3) (4) 

for field star volume emissivity jf and average pathlength A for field star radiation. The first contribution in 
the parenthesis is from external field stars, and the second is from internal field stars. The total luminosity 
of embedded young stars with volume emissivity jy is 

Ly = {4/3)7rR^jy. (5) 

The ratio Ly/Pf has to be large for the cloud temperature to increase significantly as a result of embedded 
star heating. This ratio gives a critical ratio of volume emissivities 

Z^ > 1 + ^ ^ 102 _ 103 (6) 

jf R 

for embedded star heating; here we have taken an external path length equal to 10 — 100 times the cloud 
size, considering this as the ratio of cloud mean free path to size, or the inverse of the volume filling factor 
of these ultracold CO clouds. 

The volume emissivity of embedded stars is jy = (L/M) ep for average stellar luminosity-to-mass ratio 
(L/AI) , star formation efficiency e (the ratio of the star mass to the total cloud mass), and gas density p. 
To find the average {L/M) as a function of time, we integrate the pre-main sequence stellar luminosities 
over the stellar mass function, 

L\ Q L{M,t)n{M)dM 

m)^' ^^l'^ Mn{M)dM 

using the Alexander -f- Rogers & Iglesias, CM model for £(M, t) in D'Antona & Mazzitelli (1994) and 
n{M)dM — niog{M)d\ogM from equation^. These pre-main sequence models are only for M < Mjj = 2.5 
Mq, so our results are valid only in this limit. Equation gives {L/M) for a time t after a short burst of 
star formation. To get (L/M) from continuous star formation that has lasted for a time t, we use: 

M„ 



L\ .._ /o dt J^l;; L{M, t)n{M)dM 



^^ ) y,cont /o dt /^; MTZ{M)dM 

for star formation rate by number, TZ{M), which is proportional to n{M) for a uniform star formation 
rate. The lower limits to these integrals are taken to be Mj instead of the absolute minimum stellar mass, 
because the form of the IMF is not well known below the peak and because lower mass Brown Dwarfs 
should not contribute much to the luminosity anyway. 

The burst and continuous {L/M) are shown in Figure 1 as functions of time since the burst and 
as functions of age, respectively, for Mj = 0.01 M©, and for mass ranges 0.02-1 and 0.02-2.5 M©. For 



continuous star formation, {L/M) ^ 2 L©/ M© ^ 4 erg s ^ g '^ for this IMF. For the larger mass range, 
the burst (L/M) levels off because massive stars dominate the light and reach the main sequence after 



10 years. 
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The volume emissivity of field stars is jf ~ {L/M)rpf. For Sb galaxies like M31, [L/M],- ^ 0.5 
Lq/Mq (Roberts & Haynes 1994). Thus the condition for significant young star heating is 

^ (10^ - 10^) {L/M)^ 25 - 250 

Pf "" e {L/M\ ^ e ■ ^ > 

The background stellar density in the inner M31 disk is several Mq pc~^, and the density in the ultracold 
CO clouds is in the range 1-100 Mq pc~^, given above, so p/pf ~ 0.5-50. This is comparable to the critical 
density ratio for young star heating only if e >> 0.1. 

We conclude that a Brown Dwarf + stellar population with M < 2.5 Mq would not heat the ultracold 
CO clouds in MSI noticeably if the star-to-cloud mass ratio is less than 10%. 

Another way to get a limit on how many Brown Dwarfs and other young stars might be present is from 
a limit on the mass of the most massive star that can be associated with these clouds. The largest stellar 
mass, Mmax, that is likely to come from an IMF with a total stellar mass Mtot is given by 

/- Mn{M)dM 

Figure 2 (solid lines, left axis) shows the total stellar and Brown dwarf masses versus Mmax (the two lines 
nearly overlap). If the largest pre-main sequence or stellar mass that can be hidden in and around the 
ultracold CO clouds in M31 is ~ 2 Mq, then the total mass of all the Brown Dwarfs+stars has to be less 
than ^ 50 Mq. Figure 2 (dashed lines, right axis) shows the numbers of Brown Dwarfs (upper dashed line) 
and stars as functions of Mmax- With Mj ~ 0.01 Mq, there could he ^ 10'^ Brown Dwarfs and ^ 80 stars 
less massive than 2.5 Mq in each of the large ultracold CO clouds in M31. 



6. A low star formation rate 

This number of ^ 10^ hidden Brown Dwarfs in ultracold CO clouds may seem large compared to what 
has been found so far in the Solar neighborhood, but the associated number of normal stars and the total 
mass in the form of stars or Brown Dwarfs is remarkably small for such a cloud with an estimated ^ 10® 
Mq of gas (Loinard & Allen 1998). If the efficiency of star+Brown Dwarf formation is anything like it is 
locally, namely > 1%, then so many stars should have formed that the IMF would have sampled out far 
enough into the high mass tail to produce 0-type stars, which would be seen easily. This implies that star 
formation is unusually inefficient in the M31 ultracold CO clouds. 

One possible explanation for this is that the clouds are much less dense than the excitation densities 
given by Loinard & Allen (1998). These authors found, on the basis of excited-state CO line ratios, that the 
gas is optically thick in the 3-2, 2-1, and 1-0 ^^CO lines. In that case, the density required for excitation is 
less by the factor 1/t for opacity t than the density that makes the collision rate equal to the spontaneous 
transition rate A. Thus the molecular density can be ~ 10 times less than their estimate, and the clouds 
much less strongly self-gravitating. If the clouds arc like galactic translucent clouds, they may form no stars 
at aU (e.g.. Hearty et al. 1999). 



We can use the example provided by M31 to estimate the importance of Brown Dwarf formation 
generally. Suppose the efficiency of Brown Dwarf+star formation in ultracold gas is the upper limit given 
by the M31 clouds, which is ~ 10"'* in 10^ M© clouds to avoid stars more massive than 2 Mq. Suppose 
also that half the stellar-like mass goes into Brown Dwarfs and that the ages of the ultracold clouds are 
at least ^ 10^ years, which is a modest fraction of the shear time in the inner disk M31. Then the Brown 
Dwarf formation rate is < 5 Mq My"-'^ per 10^ Mq of ultracold gas. This implies a gas consumption time 
that is very large, 2 x 10*^ years, which means hardly any conversion of gas into stars or Brown Dwarfs in 
a Hubble time. To get a significant mass in Brown Dwarfs, the efficiency of star formation has to be larger 
than this by a factor of ~ 100, but then the resulting luminosity would heat an ultracold cloud to normal 
temperatures, shutting off Brown Dwarf production in favor of regular stars. 

The formation of Brown Dwarfs in ultracold gas requires a negligible mass of normal stars so the gas 
remains cold. This would occur naturally in the lowest mass clouds from an ensemble of clouds if the stars 
formed randomly. For example, if the maximum total mass that can form in a cloud before a massive star 
is likely to heat the gas is ~ 100 Mq, from figure 2, and if the efficiency is large, 10%, then the clouds that 
make Brown Dwarfs have to be smaller than 10^ Mq. Such small clouds are rarely isolated, however, and 
nearby gas can still make massive stars that provide heat. Thus is it unlikely Brown Dwarfs can form with 
a significant total mass if the IMF is simply shifted towards a lower peak mass with a Salpeter spectrum 
above this. An upper mass cutoff in star formation, which is not present in the theory, would seem to be 
necessary if Brown Dwarfs are ever found to dominate the mass in a region. 



7. Conclusions 

The ultracold clouds in M31 provide a good test for theories of the IMF. Because of their low 
temperatures and normal-to-high pressures, these clouds should have a thermal Jeans mass of only ^ 0.01 
Mq, making the IMF shift towards lower mass with the same slope at higher mass. In that case, ultracold 
clouds should produce half of their stellar-like mass in the form of Brown Dwarfs. If the total number of 
such Dwarfs per cloud is < 10'^, then they and their accompanying H-burning stars would not significantly 
heat the cloud, nor be visible in existing surveys. A deep K-band search for Brown Dwarfs in ultracold gas 
would be necessary to see them. Considering the luminosities of such pre-main sequence stars found locally, 
which is iiT ~ 13 mag dereddened (Luhman et al. 1998), the K band luminosity of such a star in M31 would 
be ~ 21 mag plus extinction. 

Helpful comments by the referee are appreciated. 
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Fig. 1. — The luminosity to mass ratio, in units of Lq/Mq, from pre-main sequence and main-sequence 
stars in the mass ranges indicated, are plotted as functions of time for a burst at i = and for continuous 
star formation. These L/M ratios assume an initial mass function of the form given by equation (3) with 
Mj = 0.01 Mq. 
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Fig. 2. — The total stellar and Brown Dwarf masses are shown with the left hand axis as solid lines (these 
masses are approximately equal, so the two lines are nearly superposed), and the cluster stellar and Brown 
Dwarf numbers arc shown with the right-hand axis, as dashed lines. The values assume Mj — 0.01 M©. 
The Brown Dwarf dashed line is above the star dashed line. Both are plotted as functions of the mass of the 
largest star that is likely to form. 



